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3. 


Summary. 

A  method  to  calculate  the  electromagnetic  scat¬ 
tering  properties  of  a  cluster  of  spheres  of  arbitrary 
radia  and  (possibly  complex)  refractive  indexes  is  pro¬ 
posed.  The  approach  takes  proper  account  of  multiple 
scattering  effects  through  an  appropriate  addition  the¬ 
orem  for  Vector  Helmholtz  Harmonics  which  is  preliminar- 
ly  formulated.  No  approximation  is  required  but  for  the 
truncation  of  the  multipolar  expansion  of  the  scattered 
field.  Group  theory  is  also  used  to  factorize  the  resul¬ 
ting  system  of  linear  nonhomogeneous  equations. 


i 


•<* 


1 .  I  ntr oduct  i  on.. 

As  is  well  known,  the  problem  of  the  electromag¬ 
netic  scattering  from  irregularly  shaped  objects  has  not 
been  soved  in  general  but  from  a  conceptual  point  of  view^ 
For  this  reason  this  paper  dels  only  with  the  scattering 
from  spherical  objects  or  from  objects  composed  of  sphe¬ 
rical,  though  not  necessarily  homogeneous  scatterers.  We 
shall,  in  fact,  introduce  a  model  scatterer  composed  by 
a  cluster  of  spheres  whose  relative  positions,  radia  and 
(possibly  complex)  refractive  indexes  are  assumed  to  be 
known.  We  are  able  to  describe  the  field  scattered  by  such 
an  object  in  a  rather  simple  way  and  to  take  account  of 
multiple  scattering  among  the  spheres  in  the  cluster.  It 
will  become  apparent  that  such  a  model  scatterer  should 
be  suitable  to  approximate  the  scattering  properties  of 
molecules,  even  the  most  asymmetrical  ones. 

In  the  course  of  our  description  we  wi  I  I  make 
large  use  of  the  vector  solutions  of  Helmholtz  equation 
in  spherical  coor-ji  nates.  Although  these  functions  are 
thoroughly  described  in  the  literature,  the  material  is 
rather  scattered.  Therefore,  in  the  first  few  paragraphs, 
we  summarize  the  main  properties  of  these  functions.  Of 
course  our  description  does  not  pretend  to  be  complete, 
for  it  will  be  restricted  to  the  topics  of  interest  for 
our  purposes.  Further,  we  shall  give  a  detailed  descrip¬ 
tion  of  our  model  scatterer  and  discuss  the  techn i%ne-{to 
calculate  the  scattered  field  and  related  quantities. 
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2.  The  field  equations. 

Any  theory  of  the  electromagnetic  scattering 
should  start  from  Maxwell  equations  for  stationary  media, 
which  we  shall  rewrite  here  in  gaussian  units  ' 


7xH  s  ££  j  4. 

Vx  E  =  -  i-  J 1 

d  s  ili r  ^ 

7»  B  =  o 


(1-40 


U-U) 

ti-ic) 

U-4d) 


together  with  the  constitutive  equations  of  the  medium 

D  a  e  £  U-iO 

U-u>) 

j  s  (T  E 

In  the  f  o  I  lowing  we  shal  I  deal  with  nonmagnetic  l) 

and  isotropic  media  and  assume  an  harmonic  time  dependen¬ 


ce  of  all  fields  so  that,  <|.g. 

e(x) 


U-  i) 


With  these  assumptions,  equations  (2-la)  and  (2-lb)  read 

Vx-B  * -.CktfE  (4-4 a) 


7x  E 


=  i  k  B 


(MM 


where  k= i*V  c  is  the  magnitude  of  the  propagation  vector 


a-  5)  ' 


t 
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is  the  complex  refractive  index.  Furthermore,  eq.  (2-lc) 
becomes 

V  •  -  O  .  (1-6) 

Although  we  did  not  assume  that  the  medium  be  homogeneous 
so  that  both  £  and  can  depend  on  the  coordinates,  we 
need  not  bother  with  thir  frequency- dependence  thanks  to 
the  assumption  of  harmonic  t i me- dependence,  eq.(2-3),  pro¬ 
vided  that  the  values  of  £  and  be  those  appropriate  to 
the  frequency  at  hand. 

It  is  now  a  simple  matter  to  show  that  and  B 
are  the  solutions  of  the  equations 

V* V x  E  -  tAV  £  =0 

/>**  ^ 

v*vx  B  -  k Vis  -  .ik  vV'xE  (l-n) 

which  are  not  decoupled  on  account  of  the  assumed  non-ho¬ 
mogeneity  of  the  medium. 


3.  The  vector  Helmholtz  equation  and  its 
so  I ut ions. 

In  order  to  discuss  the  structure  of  the  solu¬ 
tions  of  equations  (2-7)  let  us  consider  first  a  homoge¬ 
neous  medium.  The  above  Equations  then  become 


(v*  +  Kl)  £  =  o  i 

(V*  +  K1)  'B  r-  o  ' 


with  K  ~  kn,  i.e.  both  t:  and  B^  should  satisfy  a  vector 
Helmholtz  equation.  This  statement  means  that  each  of  the 
components  of  E  and  B  is  the  solution  of  a  corresponding 
scalar  Helmholtz  equation.  However  this  is  in  no  way  true 
in  any  other  system  of  coordinates:  even  in  orthogonal 
coordinates  it  is  impossible  to  separate  eqs.  (3-l)  into 
three  scalar  equations  each  involving  only  one  component 
of  the  field  .  It  si  therefore  convenient  to  search  for 
general  vector  solutions  of  the  Helmholtz  equation 
a.  Hansen's  vectors. 

Let  us  consider  the  scalar  Helmholtz  equation 

=  0  (3-3) 

and  let  C  be  a  vector  operator  capab I e  of  acting  on  its 
solutions,^.  Then  we  have 

<3(V*+  K*)^>  =  (Vl+  K4)©  f  +  [S,  O-  2.) 

i.e.  for  any  vector  operator  such  that  of  is 

a  vector  solution  of  Helmholtz  equation^.  Now  two  opera¬ 
tors  with  this  property  are  the  gradient  operator^V  , 

which  is  proportional  to  the  linear  momentum  operator  P'-.t  V 

*  , 

and  the-  angu  I  ar  momentum  operator  =  •  We  have 

then  to  search  for  a  third  operator  for,  on  general  mathe¬ 
matical  grounds,  Helmholtz  equation  should  have  three 
linearly  independent  vector  solutions  for  each  value  of 
K.  1  o  this  end  let  us  recal  I  that  if  /\  is  a  vector  solu¬ 
tion,  also  VxA  is.  Thus  the  only  other'  independent  sol  u- 


8. 


t  i  on  i  s  7x7ji^  •  Ultimately,  given  any  scalar  solution  of 

“  .5,7) 

Helmholtz  equation,  the  three  vector  functions  ' 

L  =  ,  Ncivxtf  (3-4) 

known  in  the  literature  as  Hansen's  vectors  form  a  comple¬ 
te  set  of  vector  solutions.  The  factor  l/K  in  the  defini¬ 
tion  of  N  has  been  introduced  so  that 


V  x  M  =  N 


U-5) 


b.  Irreducible  spherical  tensors. 

Hansen's  vectors  are  defined  in  any  system  of 
coordinates  and  are  therefore  quite  general.  They  do  not 
form,  however,  an  ortogonal  set  for 


L •  M  =0  ,  M  •  W  =0 

—  —  I  r~  ~ 


but 


k'%  4° 

and  this  may  be  trublesome  in  several  applications.  We 
shall  therefore  introduce  a  new  independent  set  of  solu¬ 
tions  which,  although  defined  only  in  spherical  coordina¬ 
tes,  form  a  complete  orthogonal  set  and  have  very  useful 
mathematical  properties.  To  this  end  let  us  search  for 
the  tr ansf ormat i on  properties  of  a  vector  field,  £(x,y,z), 
under  infinitesimal  rotations  of  the  coordinate  axes.  As 
is  well  known,  the  infinitesimal  rotation  operator  is' 

A 

"R(«i«)  =  (4  +  *  4%  •  J  ) 


.  8) 


A  * 

where  J_  is  the  angular  momentum  operator,  so  that  *\:  t 

u  +  *  (1-6) 


s 
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Now,  by  working  the  tr ansf ornrnt i on  of  the  rect jngul or  com¬ 
ponents  of  F  and  comparing  with  eq.  (3-6),  one  easi ly  sees 
that 

A  A  / 


j=  l  +s  =  -  irxv  -t s 

x.  ~  ~ 


where  S  is  a  set  of  three  3x3  matrices  which  represent 

9) 

the  intrinsic  spin  of  the  vector  field  .  The  simultane- 

*  2  A 

ous  eigenvectors  of  S  and  S  can  be  shown  to  be 

l±l  *  +  (W) 

where  e  ,  e  and  e  are  unit  vectors  a  I  ono  the  x,  y  and 
—x  — y  ~z 

z  axes,  respectively.  It  is  apparent  from  the  above  equa¬ 
tions  that  the  intrinsic  spin  of  any  vector  field  is  1. 

2 

As  the  simultaneous  eigenvectors  of  L  and  L  in  spheri¬ 
cal  coordinates  are  the  well  known  spherical  harmonics,  Y^, 
we  can  couple  the  Y's  and  the  ^  's  by  means  of  the  vector 
coupling  coefficients  and  get  the  simultaneous  eigenvectors 


of  J  ,  J  , 


L2  and  S2 


Tji/f)  =  (i-8) 

M 

The  vectors  Xj^(r)  ore  the  components  of  an  irreducible 
spherical  tensor  pf  rank  2J+1  and,  under  rotation  of  the 
coordinate  axes  transform  according  to^^ 


i 


(3-9) 


while  their  parity  is  (-)  .  The  triangular  condition  on 
the  C I cbsch-Gordan  coefficients,  imposes  that 


10. 


for  each  value  of  J  there'  exist  only  three  mutually  ortho¬ 
gonal  tensors,  ^  .  Now  it  i  s  of  immediate  ver  i  f 
cation  that  the  vector  functions 


Aju(ri  =  (3-1°) 


where  f  i s  a  spherical  Bessel,  Neumann  or  Hankel  func¬ 
tion,  are  solutions  of  the  Helmholtz  equation  in  spheri¬ 
cal  coordinates.  They  will  be  referred  to  in  the  fol lo¬ 
wing  as  Vector  Helmholtz  Harmonics  (VHH). 

c.  Vector  Spherical  Harmonics  (VSH). 

When  one  has  to  impose  to  the  e I ectromagnet i c 
field  the  customary  boundary  conditions  at  the  surface 

of  a  sphere,  the  VHH's  are  not  the  most  useful  set  of 

M 

basis  functions  for  T  are  neither  tangent  nor  ortho- 

gonal  to  the  surface,  whereas  T  is.  We  can,  however, 

M  ^ 

take  linear  combinations  of  Zll+1  w^*  • d°  are  orthogo¬ 
nal  and  tangent  to  the  surface  of  the  unit  sphere.  There¬ 
for*  we  def i ne 


&«<*>  =  - 

■  Z!Ur) 

(l-Mo,) 

i\„(:)  =  - 

Jjzii.  + 

|U+i  -t-L  +  l  * 

llsrrr 

•T*  M 

-Ill-1 

(VdU) 

rr 

th 

A  L  L  -<i 

The  above  functions,  known  in  the  literature  as  Vcdt^o** 
Spherical  Harmonics,  arc  orthogonal  to  each  other  andj 


11. 


when  multiplied  by  f^(kr)i(  form  a  complete  set  of  solu¬ 
tions  of  the  vector  Helmholtz  equation.  Thir  relation 
to  Hansen's  M  and  N_  vectors  in  spherical  coordinates  is 
given  by  the  equations 


We  notice  that  the  solenoidal  character  of  M.  ..  and  of 

— LM 

— -LM  '  's  £lui'te  evident  from  the  above  equations, 

2 

allows  to  expand  any  solenoidal  field,  such  as  n  £  and 
B,  in  their  terms  only.  In  this  case  ^LM  '  of  par i ty 
(-)*”,  is  said  to  represent  a  magnetic  multipole  of  order 
2^,  while  ,  of  parity  (-)^  represents  an  electric 

multipole  of  the  same  order.  The  relation  between  Hansen's 
J.  vector  and  the  VSH's  is  rather  complicated  and  will, 
however,  never  be  used  in  our  work. 


4.  An  addition  theorem  for  VHH's. 
a.  The  addition  theorem. 

In  the  course  of  our  work  we  shal I  need  to  re¬ 
late  to  each  other  the  VHH's  centered  at  the  eigin  of  two 
mutually  translated  systems  of  spherical  coordinates.  To 
do  this  we  start  from  the  addition  theorem  for  scalar  Helm 

holtz  harmonics  wch i ch  we  rewrite  here  in  a  form  slightly 

12) 

different  from  that  reported  by  Nozawa  : 


12. 


where  the  quantities 

with  r#  =  r  -  Rf  are  the  matrix  elements,  in  the  angu I ar 
momentum  representat i on  of  the  free  space  propagator  for 
spherical  waves.  In  eq.  (4-1)  and  (4-2)  when  f^=  ,^=j^ 

and  g^,=  Jl>  •  but  when  f^  =  h^ 

+>=  Jx  -  ■  ’r'>* 

and  the  quantities 

I^i'/CjLH)  =Jy*n.x„  Yx„.„m  JA  «-3) 

13) 

are  the  well  known  Gaunt  integrals  .  Now  we  recall  that 
a  VHH  is  defined  as 

AjVn-fjWT/lcru 

=  ZC(-l,i.,J;y  ,M+h)|1_(kr)YLfl^,(T:)  (4  -4^ 

so  that  we  can  substitute  eq.  (4-l)into  eq.  (4-4)  to  get 

Ajl (jc )  3  X- J ) X.  i. h+/* )  * 

(?')  fyc 

which  can  be  written  as 


**♦***,7'.'’  : 


A^<t)  =  Z C« «.  J i -/ . "vOl  &lVl 


ZjCU.l’.j';  ysM")  t",J“ (£') 

3* 


through  the  use  of  the  inverse  to  eq.  (3-8)*  If  we  now 
put  M '  =  W'-J*  and 

we  get 

14) 

which  is  the  required  addition  theorem 

Equation  (4-5)  can  be  specialized  to  Hansen's 
M  and  N  vectors  through  the  use  of  the  relations**^ 

nLh  =  {Win  =  *[$SiL<TLML.t -jiZu  («1"J  ■ 


Indeed  we  have  for 


(i-i) 


I 
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whence,  on  account  of  tho  d i vergence I  ess  character  of 
the  recursion  relation  follows 

®  C.u  W  (‘-?) 


which  can  also  be  proved  by  direct  calculation  making  use 
of  the  recursion  properties  of  the  C I ebsch-Gor dan  coeffi¬ 
cients.  With  the  help  of  equation  (4-7),  eq.  (4-6)  can  be 
put  into  the  f orm* 

MLnM-  L  ^(f)  *  M-») 

l'm*  l 

where  we  put  ______ 

I  a  .  JilLd  e*'"  «- 5) 

~Jl'l',LL  1  llh'LH  k  0L'L+i,Ll 

The  functions  M  and  are  identical  to  M  and  N,  respecti¬ 
vely,  but  for  the  substitution  of  g^  to  f^.  In  the  fol lo¬ 
wing  the  quantit'es  U /  and  w  will  be  indicated  by  h,LH 
and  .respectively,  when  contains  h^  and  by  ‘in 

and  o^l'hUh  .respectively  when  ^  contains  . 

b.  Matrix  elements  of  the  dyadic  Green 
funct i on . 

In  this  section  we  will  show  that  the  quantities 
arc  the  matriw  elements  of  the  dyadic  Green  function 
for  free  space  propagation  of  spherical  vector  waves.  To 
this  end  we  recall  that  this  function  is  the  solution  of 
the  inhomogeneous  vector  Helmholtz  equation 

(?l+  kl)G(r,r’)  =  -  4-rr4S(r-r')  , 


15. 


and  has  the  form 


G0C)r')  =  l 


gl  k  I  £- 

lr-  r'l 


(4.40) 


where  is  the  unit  dyadic.  If^for  greater  generality  we 
consider  two  points  in  space  at  and  Rp  and  put  ; 

V,-yR^  »cq.  (4-10)  can  be  rewritten  as 

G(r«,rjU-4  U.  U) 

with  .  Neuman  expansion  of  eq.  (4-ll)  is 

fit*,!*)  .Ml  aLU^-^Xn(X<-ty IkNlXfojt) 

on  the  assumption  that  .  By  assuming  further 

that  /iK'fUft  the  addition  theorem  of  eq.(4-l)  can  be  ap- 
pl  ied  to  •<>  t  ( 1 1 f<<  ‘  '£*«<a0  YtM  (T<  *£«/0  * 

where  wc  expanded  the  unit  dyadic  in  a  spherical  basis: 

Now  we  recall  that  the  spherical  harmon i cs  and  the  irredu¬ 
cible  spherical  tensors  arc  related  through  the  equation 

\T  Xh  (r ).  Z  c  a,  1 ,  J;  r,H)  r'^Ct ) 

3" 

so  that  eq.  (4-12)  can  be  rwritten  as  { 

~ (U,Tt,)  '  /,T  ;  TjT'^tJO  C(l,LJ;y  ,11)  • 


U-UI 


16. 


-/, «')  t.U) 


which,  through  the  position  M  - J*  =  m,  M'-yit=m',  takes  the 
final  form 

S<*  W  =£££  j^>  Ti^>  (^5) 

Eq.  (4-13)  shows  that  the  quantities 

®*  A  (4“  I 

are  just  the  matrix  elements  of  G  with  respect  to  VHH's 

centered  at  different  sites  of  space. 

We  want  to  remark  that  the  same  result  would  be  obtained 

if  we  assumed  1  if*  o^j ) $■  rp  and/or  but  for  that 

the  VHH's  centered  at  JR  ^  would  contain  h^  instead  of 

and  G.  , ......  would  contain  j  instead  of  h.  .  Finally  we 

L  M  LM  >  A 

notice  that  the  quant i t i esdef i ned  in  eq.  (4-14)  are  the 
off  diagonal(in  the  site  indexes)  elements  of  G,  whereas 
the  quantities  defined  in  eq.  (4-5)  are  the  corr espond i ng 
on  diagonal  elements. 


5.  Scattering  from  a  radially  symmetric 
sphere,  •, 

The  theory  of  electromagnetic  scattering  from 

a  homogeneous  sphere  is  well  known  since  the  work  ?if 

1/)  '7  - 
Mic  .  I  n  this  section  we  wi  I  I  study  the  scattering 


from  a  radially  symmetric  sphere,  i.e.  with  n=n(r),  of 
radius  b,  with  the  twofold  purpose  of  finding  the  general 
expression  of  E  and  B  within  an  inhomogeneous  medium  and 
of  defining  several  quantities  which  wi  I  I  be  useful  for 
subsequent  work. 

a.  The  electromagnetic  field  within  a 
radially  symmetric  sphere 
As  shown  in  section  2  the  E  and  B  fields  within 
a  nonhomogeneous  medium  should  satisfy  eqs.  (2-7).  However 
in  such  a  medium  V*  Sf  «  so  that  eq.  (2-7a)  cannot  be 
cast  in  Helmholtz  form.  It  is  therefore  more  convenient 
to  search  for  solutions  of  the  Maxwell  equations.  To  this 
end  we  expand  I:  and  B  in  a  series  of  VSH's 


LM 


Lh 


(5-U) 

(5-lM 

which  satisfy  eqs.  (2-6)  and  (2-ld)  ,  respectively,  for 
any  choice  of  the  radial  functions  R^  and  S^.  These  latter 
arc  determined  by  substitution  of  eqs.  (5-l)  into  the  Max¬ 
well  equations  (2-4) «  Indeed,  it  is  easily  seen  that  the 
radial  functions  should  be  the  solutions,  regular  at  r~0; 
of  the  equations  '' 

f  j£  -  illtlL  +  kVl(T«t(^))=  o 

**  %-r 

[A  -  ilillL  jeL  -v  kV| (TSt(o)  =  o  ,  M*>) 

*-4rl  \x  dr  J 
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respectively.  The  expansion  coefficients  C.  and  D.  are 

LM  LM 

then  determined  by  the  requirement  that  £  and  B  undergo 
the  appropriate  boundary  conditions  at  the  surface  of  the 


med i um. 


b.  The  incident  plane  wave. 


We  assume  that  the  field  inciding  on  the  sphe¬ 
re  has  the  form  of  a  circularly  polarized  plane  wave  of 


wavevector  k: 


-  ,  .  „  i  Lr 

.1 

ifi.  =  Ai  (fi,  o  * 


*5^  =  fail) 

where  e^  and  are  unit  vectors  orthogonal  to  k  and  to 
each  other  and  according  to  the  polarization.  The 

direction  of  the  propagation  vector,  k,  is  assumed  to 
be  quite  general.  Now,  if  we  assume  for  JE^  and  the 

expansions  in  VSH's 

-L  [a^LM  jLMX LH(i)  +  ^  7xjt(kOXLM(t)] 

l  Qy  Z_[  ty M  jt(W  XLH(? )  +  a^LH  4j- 7xjt(U)  Xuh(jf )] 

3) 

a  strai ghtf orword  but  lenghty  calculation  shows  that 

aAjLH  *  + 

-i  a  ^ 

Thus  the  expansion  of  the  incident  wave  in  VSH's  is 

=Z.^tft(k)[j1.(kr)XLM(i)  +  ^^  VxjL(kr)Xt#i(t)]  (5-Ao.) 
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Vfn  *!  <«U 


where 


W*jLM(ic )  =  h riLt&i  +  *  ^g*)‘ 


(5*2) 


contains  the  dependence  on  the  direction  of  k. 

Cm  The  scattered  field. 

The  only  ingredient  we  need  to  solve  our  scat¬ 
tering  problem  is  the  scattered  field  which  we  assume  in 
the  form 

=  £_  [-By l*(£)  +a^h i  v  Xl„(X)]  (5-^1 


where  the  Hankel  functions  of  the  first  kind,  h^(kr),  en¬ 
sure  the  correct  behavior  of  the  field  at  infinity:  eqs. 
(5~6),  indeed,  contain  only  outgoing  spherical  waves. 

Now  we  rccal I  that,  at  the  surface  of  the  sphere,  the 
tangential  components  of  £  and  B  should  be  continuous 
as  well  os  the  radial  component  of  B,  while  the  radial 
component  of  £  should  satisfy 


A 

r 


where  the  index  1  refers  to  the  interior  of  the  sphere. 


By  imposing  the  above  boundary  conditions  we  get,  fcr 

each  L,M,  six  equations  among  which  C.  and  D.  ..  ,  the 

LM  LM 

coefficients  of  the  internal  field,  are  easily  elimina¬ 


ted.  The  remaining  equations  then  yield  for  A and 

B>jLM  ocluo^'ons 


=  W, 


. «tr 


r-  o 


yin 


^  (rAL(kr)J  JT=  b 


=  1WY«^- 

Equations  ( 5—  7 )  solve  completely  our  problem  since  all 
the  quantities  of  interest  can  be  expressed  in  terms  of 
the  A  and  B  coefficients.  In  particular  the  scattering, 
absorption  and  total  cross  sections  are  given  by 

Ay„|1  +  (*-*«•) 

r,.U)  £[*_  |AV(1  +  ,{|*-  I6>|in+l|‘] 

»  k  Lh 

*')'  =  *^T  L.  SU  (^/VJUH  +  (5-Sc) 

where  fit  denotes  the  real  part  .  Note  that  for  a  sphere 


I 

4 


the  dependence  of'  the  A's,  B's  and  CTrA  on  the  polarize- 


t  i  on  index,  ,>j  .  However,  we  will  see  in  the  next  section 
that  the  cross  section  of  a  cluster  of  spheres  can  st i  I  I 
be  cast  in  the  form  of  eqs.  (5-8),  but  with  on  actual  de¬ 
pendence  on  the  polarization  of  the  incident  wave. 


6.  Scattering  from  a  cluster  of  spheres. 

As  is  we  I  I  known,  the  scatterers  in  the  most 
common  aerosols  are  far  from  spherical  and  their  proper¬ 
ties  cannot  be  described  by  Mie  theory  but  to  low  appro¬ 
ximation.  We  shall,  therefore,  introduce  a  model  scatte- 
rer  whose  features  allow  to  fit  the  properties  even  of 
highly  asymmetric  objects,  although  their  scattered  field 
can  be  calculated  without  too  much  computational  effort. 

We  define,  in  fact,  our  model  scatterer  as  a  cluster  of 
N  nonmagnetic  spheres  whose  centres  lie  at  R ^  and  whose 
radio  and  (possibly  complex  and  radially  symmetric)  re¬ 
fractive  indexes  are  b^  and  n^  ,  respectively.  The  field 
inciding  on  the  cluster  is  assumed  to  be  the  plane  wave 
described  by  eqs.  (5-4)/  wh i  le  the  scattered  field  wi  I  I 
be  written  as  _ 

6';’  =iL[A*t in 

* =1 L  [<„  A“lKi-vx!L(U,)xlB<r,)]  (c- U) 

I  IH 

with  i.e.  as  a  superposition  of  the  fields 

0( 

scattered  frm  the  single  spheres.  The  coefficients 
and  ,  however,  are  not  merely  given  by  eqs.  (5-7) 
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but  ore  calculated  so  as  to  account  for  multiple  scatte¬ 
rings  among  the  spheres.  In  other  words,  the  field  inci¬ 
dent  on  theflf-th  sphere  is  assumed  to  be  that  of  the  in¬ 
coming  plane  wave  plus  that  previously  scattered  by  all 
other  spheres.  This  effect  is  achieved  by  rewriting  both 
the  incident  and  the  scattered  fields  in  terms  of  VSH's 
centered  at  through  the  addition  theorem  of  sect.  4- 
We  get 

V<jt,UO  A)j 

+  4j^[<^LiM,lh  ^  Vx  &„(*>]} 

£^=2l{^LM  X1>H(5»)+'&/JU,  ^VxiL(^)Xt*  &*) + 

+^tn6K^,thjJk)X,.v&)  VxjcdrJ  &!».  (&))]] /‘'V 


and  analogous  equations  f  or  a’B!  and  i  G>*  \  The  field  within 

*1  ~1 

the  pC-th  sphere  is  assumed  in  the  form  of  cqs.  (5-6) 
with  r^  substituted  to  r. 

Now,  wc  take  the  dot  product  of  cqs.  (6-2),  (6-3)  and 
(5-6)  in  turn  w  i  th  ) ,  tfjJLU  )  and  It*  X^jL  )  *"d 

get  the  radial  and  tangential  components  of  the  fields 
at  the  surface  of  the  d-th  sphere.  Imposition  of  the 
boundary  conditions  and  integration  over'  the  angles  yield 


23. 


for  each  Of,  I  ,  m,  six  equations  among  which  Ch  and 

*\in ) 

the  coefficients  of  the  in ter no!  field  ore  cosily  elimi¬ 
nated.  This  possibility  allows  to  get,  for  each  d,\,m, 
two  equations  involving  only  the  A's  and  B's  as  unknowns 

=  "^5  ^uh  (i)  «-<*) 

£fJ<Ve  V  -  *&.<■!  = 

=  -Z.WyB(C)Q^lB  U-H) 

In  eqs.  (6-4)  we  define 

+/>)  3/1 IH 

while  ($^  and  are  st  i 

with  obvious  change  of  the  arguments.  The  system  compo¬ 
sed  of  equations  (6-4)  for  all  values  of  fltl,m,  allows 
a  complete  and  uni  cue  determi not i on  of  the  scattered 
field  and  thus  of  all  the  quantities  of  interest.  As 
an  example,  to  calculate  the  cross  sections  of  the  clus¬ 
ter  we  rewrite  the  scattered  field  in  terms  of  VSH's 
centered  at  a  single  site, say  R 0 ,  through  application  of 
the  addition  theorem  of  sect.  4.  We  get  for  the  electric 
field 


defined  as  in  eqs.  (5-7) 
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\  VxijdiOX,.,,, 


(4-6) 


with  r0  —  £  -  R0  ,  and  an  analogous  expression 
Now,  by  defining 


ZV«-LL[^  C. 
-B  -  7  7  iv 


for  • 


the  scattered  field  can  be  written  in  the  form  of  that 
scattered  by  a  single  sphere,  so  that  the  cross  sections 
of  the  cluster  are  still  given  by  eqs.  (5-8)  with  A's 
and  B's  substituted  for  A's  and  S's,  respectively. 

It  is  to  be  noted  that,  since  the  cluster  leeks  the  sphe¬ 
rical  symmetry,  the  cross  sections  are  actually  dependent 
on  the  polarization  index, Aj  .  Furthermore  we  notice  that 
the  choice  of  the  point  R0  is  quite  arbitrary  but  that 
according  to  sect.  4,  eq.  (6-6)  is  valid  for 

**  '  *v  "  * 

for  any  oL  ,  i.e.  in  the  region  external  to  a  sphere  cen¬ 
tered  at  Rc  and  including  the  whole  cluster.  The  radius 
of  this  sphere  can  be  minimized  by  choosing  R0  at  the 
centre  of  mass  of  the  cluster,  but  the  choice  is  not  cry- 
tical  for  the  scattered  field  is  always  observed  at  dis¬ 
tances  I arger  than  thg  size  of  the  cluster. 


* 
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7.  Con vergency . 

The  method  described  in  the  preceding  sections 
does  not  require  any  approximation  but  for  the  mu  1 1 i po- 
I ar  expansion  of  the  scattered  field,  eqs.  (6-1).  It  is 
therefore  of  fundamental  importance  for  the  feasibility 
of  the  method  itself  to  discuss  the  rate  of  convergence 
of  these  expansions.  To  this  end  let  us  rewrite  the  sys¬ 
tem  of  eqs.  (6-4)  matrixwise  as 


<St-‘+  H  X 

A 

f 

ft 

Q 

(M) 


Equation  ( 7— 1 )  allows  to  identify  the  matrix  on  its  left 

hand  side  as  the  inverse  of  the  e I ectr omagnet i c  T-matrix 

2) 

of  the  whole  cluster  .  Of  course  it  is  not  diagonal  on 
account  of  the  lack  of  spherical  symmetry  of  the  cluster 
as  a  whole.  However,  the  matrices  and  ^ 

have  an  interesting  physical  meaning  of  th i r  own  and  their 
anal i cal  behavior  determines  the  rate  of  convergency  of 
the  scattered  field,  as  we  are  going  to  discuss. 

The  diagonal  matrices  (R  and  J  are  the  direct  sum  of  the 

-v  r* 

matr i ces  <$.  and  /  defined  according  to  eqs.  (5_7), 

/■V  /V 

which  account  for  the  scattering  power  of  the  C(  -  th  sphe¬ 
re  in  the  absence  of  any  other  scatterer.  Thu* the  presen¬ 
ce  in  the  cluster  of  more  than  one  scatterer  is  accounted 
for  not  oil  I  y  by  the  and  withyfl^d  but  also  by  the 

matrices  ${  and  5^  which  couple  to  each  other  all  the 


I 


spheres  in  the  cluster.  Indeed,  os  shown  in  sect.  4 ,  M 

an  ore  the  matrix  elements  of  the  free  space  dy¬ 

adic  Green  function  ond  thus  describe  the  propagation  to 


sitetf  of  the  waves  scattered  by  site  ^  .  This  remark  ful¬ 
ly  justifies  our  previous  statement  that  all  multiple 
scattering  processes  are  accounted  for  in  the  present 
theory . 

As  regards  the  rate  of  con vergenney ,  we  remark 
that  it  is  expected  to  be  fairly  good  even  when  the  clus¬ 
ter  is  not  small  in  comparison  with  the  incident  wave- 
lenght,  provided  kb^«  1  for  any  o(  .  Under  this  condition, 

and  decrease  rapidly  with  increasing  /  so  that 

«  * 

and  di  are  quite  sufficient  to  describe  the  field 

scattered  by  the  oL  -th  sphere  in  the  absence  of  any  other 

17,19) 

scatterer,  even  when  n ^  is  not  close  to  unity  ' 

The  rate  of  convergence  thus  depends  on  the  behaviour 

and  .  Now,  one  easi ly  sees  from  their  de¬ 

finition,  that  th i r  order  of  magnitude  is  determined  by 
the  Gaunt  integrals,  ,  and  on  the  spherical  Hankcl 

functions  h ^(kR^^) .  The  -integrals  do  not  vanish  only 

fdr  lL-{|5.X$  L+  l  and  decrease  very  rapidly  with  in- 
creasing  \  •  Thus,  although  the  imaginary  part  of 

tends  to  increase  when  \  the  even¬ 

tual  effect  is  to  decrease  the  magnitude  of  both 

wA 

and  with  increasing  /  ,  L  anil  This  behaviour 

is  to  be  expected  for,  when  the  intersphere  di  stanoe  Jfs 
very  largo,  the  present  theory  should  reduce  to  that  of 
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the  scattering  from  N  spherical  scattercrs  without  any 
multiple  scattering  effect.  As  a  consequence  it  is  rea¬ 
sonable  to  expect  that  our  approach  converge  we  I  I  by 
truncating  the  expans i on,  eqs.  (6-1),  of  the  scattered 

field  at  L=L  =3.  However,  the  order  of  the  system  (7-1) 

M  2 

is  d,.=  2N(L..+l)  -2N,  so  that  we  have  d  =30N,  a  very  high 
MM  3 

number  even  for  small  clusters.  Nevertheless,  if  the 
cluster  posses  symmetry  properties,  as  is  the  case  for 
actual  molecules,  we  shall  use  group  theory  to  put  the 
system  (7-1 )  in  factorized  form,  as  will  be  shown  in  the 
next  sections. 


■ 

k 


8.  Symmetr i zat f on . 

We  assume  that  he  cluster  is  left  invariant 
by  the  transformations  of  a  group  ^  of  order  g.  The  ef¬ 
fect  of  these  tr ansformat i ons  is  to  permute  among  them¬ 
selves  the  spheres  in  the  cluster  but,  in  general,  not 


all  the  spheres  are  linked  to  each  other  by  an  operation 
of  ^  .  Therefore  we  partition  the  cluster  into  sets 

of  spheres  which  are  transformed  to  each  other  by  the 


operations  of  the  group.  Of  course  this  does  not  imply 
a  renumbering  of  the  spheres  but  only  that  we  associate 
to  each  site  index,  of  ,  the  appropriate  set  index,  5” . 

In  order  to  get  the  system  of  eqs.  (6-4)  in 
factorized  form  we  have  to  expand  both  the  incident  and 
the  scattered  field  in  terms  of  comb i nat i ons  of  VSH's 
transformi ng  according  to  the  rows  of  the  irreducible 


represent  at i on  of  ^  .  To  this  end  we  recoil  that,  os 

shown  in  sect.  3,  the  VSH's  transform  according  to  the 
represent  at i on  D  of  the  Pul  I  rotation  group.  There¬ 
fore,  if  S  is  on  operation  of  such  that 

SR*  -  Rp 

with  <x  and  A  in  the  same  set,  of  course,  and  0  is  the 

22)  b 

associated  operator,  then 

L  u-i*) 

provided  S  is  a  proper  operation.  When  S  i s  an  i mproper 

operation  we  must  take  account  that  f  X  and  7X  f.  X. 

L— LM  L~ LM 

have  opposite  parity,  so  that  eq.  (8-la)  must  be  multi- 

L  L+l 

pl  i ed  by  (-)  and  eq.  (8-lb)  by  (-)  and  the  argument 


of  D 


must  be  understood  as  the  proper  rotational  part 


of  S.  As  a  consequence,  to  get  the  symmctri;cd  combina¬ 
tions  of  VSH's  we  must  apply  the  projection  operators 


(8-i) 


both  to  the  magnetic  and  to  the  electric  2  -poles.  There- 
f or e  we  write 

,  *(r  fl«H 

and 

K*f,r=  L  Z.  Cl  {  VxiijUOX^tt,)  (j-iP 

for  the  combi  not  ions  of  magnetic  and  electric  niultipoles 
centered  at  sites  of  the  d’-tli  set.  The  superscripts  )>,p 
indicate  that  the  combination  belongs  to  the  p-th  row 
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of  the  \>-th  irreducible  represent  at  i  on  and  the  index  N 
recalls,  when  appropriate,  that  one  can  get  more  then 
one  set  of  basis  functions  for  a  given  L.  Now  the  field 
scattered  by  the  cluster  can  be  written  as 


with  a  similar  expression  fori  .  We  notice  that 
work  with  unitary  irreducible  representations,  as  shown 
by  the  structure  of  the  projection  operators,  eq.  (8-2), 
and  the  coefficients  GlJt  and  have  the  property 

LL(4:M)Vt^Sw ,  (8-3) 

We  have  now  to  decompose  the  incident  field  into  parts 

belonging  to  the  rows  of  the  irreducible  representations 

23) 

of  S  .  To  do  this  we  use  the  completeness  property 

and  write 


we 


where 

4  =(W3)1  4^M^’(S)jt(kv)XlK,(r)a 

»>  ^1  * 
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The  analogous  equation  forifi*!*  is  easily  obtained  through 

**  l 

the  re  I  at i on 

i  ^  E«> 

The  field  within  the  spheres  need  not  be  symmetrized  for 
reasons  that  will  become  apparent  later.  The  internal 
field  is  therefore  still  given  by 

+  c\lh{  Vx'R.C’-.O*  ,.„(£,)]  <«-«m 

At  this  stage  we  can  write  the  system  for  the  coeffici- 

i*Yr  si*l€ 

ents  V*  m,  and  lON(L  through  the  same  procedure  we  used 
in  sect.  6.  We  rewrite  eq.  (8-4)  and  (8-6)  in  terms  of 
VSH's  centered  at  the  site  ot  belonging  to  the  y-th  set 
by  means  of  the  addition  theorem  of  sect.  4-  Then  we  dot 
multiply  the  resulting  equations,  as  well  as  equations 
(8-8),  in  turn  by  £4  )  '  %tJ  t*)  and  r^x  ?*.&>  and 

get  the  radial  and  the  tangential  components  of  the  fields. 
Imposition  of  the  boundary  conditions  and  integration 
over  the  angle  yield,  for  each  >?,p,o<  ,  six  equations 

g|  J 

among  which  C*.,  _  and  D*  „  ,  the  coefficients  of  the  i  n- 

ternal  fiold^can  be  easily  eliminated.  This  circumstance 
on  one  hand  clarifies  the  i nessent i a  I i ty  of  the  symmetri- 
zation  of  the  internal  field,  on  the  other  hand  yie  I  d{>i 
for  each  V, p,  d  ,  two  equations  involving  only  the  and 

j  ^ 

the  $6  as  unknowns: 


* 
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z  £  £  (I  I  ( o  v  y  <r?  r*  a£w )  4L  *4$  * 

III  {££  ,(S* ^SyG/ftt 3&)  k&  0$  * 


♦?&.*&■  <&  <  ■ 


°  £  V1  ( [°t w  w  4 1  d*~w  j£  LJ  (MM 

It  may  see.n  that  the  system  of  equations  (8-9)  solve  our 
problem,  but,  apart  from  its  rather  asymmetrical  form, 
it  is  still  not  completely  factorized.  However  we  remark 
that  the  spheres  within  a  given  set  are  identical  to  each 
other  so  that  $*[_  and  tf\.  do  not  actually  depend  on  the 
site  indeXjO^but  rather  by  the  set  index,iT.  Therefore  we 
can  multiply  eq.  (8-9a)  by  and  eq.  (8-9b) 

and  sum  over  the  rt's  belonging  to  tjr  and  over  m.  The  or¬ 
thogonality  relations,  eqs.  (8-5),  then  yield  for  each 
^»P*T  *  the  equations 

£ [s„  [«y-\  JO-i] -  Z7  JCko fly  =  -•2y  <«-'•*> 
Z.[SrvtW  +  JO'^3  +  £  =  -  Qy/ 
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where,  for  the  soke  of  scmplicity,  we  put  rs(f, I ,n), 
ss(r,L,N),  r  '3  (s’,  I  ,  n  ' ) ,  s  '3  (f,  L ,  N  '  )  and  define 


(8*410,) 

(8- Mb) 


with  an  obvious  meaning  of  the  parameters  e,m.  The  quan¬ 
tities  and  are  identical  to 

and  ,  respect  i  ve  I  y;  but  for  the  mutual  exchange  of 

the  a's  with  the  b's.  Moreover 

(  LH  flt*.  H'  1 

Q't°LLL^Md)ar!;l)‘[  <$«&.*■  . 

I  LH 


We  remark  that  on  the  left  hand  side  of  eqs.  (B-10)  and 
(8-11)  the  superscript  p  on  and  is  missing. 

As  wi  I  I  be  shown  later,  these  quantities  are  actual  ly 
independent  of  the  row  index,  p,  which  has  accordingly 
been  dropped.  Anyway  the  systems  of  equations  (8-10) 
solve  completely  our  scattering  problem.  Indeed,  direct 
comparison  of  eqs.  (6-1 )  with  eqs.  (8-4)  yields 

/>V«  ■  '  ^Y" 


Threforc  all  the  properties  of  interest  can  be  expressed 
in  terms  of  and  ,  the  coefficients  of 

the  symmetrized  expansion  of  the  scattered  field. 


(«-//**.) 
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9.  Discussion. 

The  consequences  of  the  factorization  effected 
by  group  theory  can  be  better  discussed  if  we  write  the 
system  with  a  given  i>,p  in  matrix  form: 


r 

2>v 

s'1. 

«•  mm 

a* 

O -i) 


As  equation  (9-l)  has  the  same  overall  structure  of  eq. 
(7-1)*  the  matrix  on  its  left  is  still  the  inverse  of  the 
electromagnetic  T-matrix  for  the  whole  cluster,  or  better 
that  part  of  the  inverse  T-matrix  which  belongs  to  the 
V-th  irreducible  represent at i on  of ^ .  This  means  that, 
although  the  lack  of  spherical  symmetry  does  not  allow 
to  get  a  diagonal  T-matrix*  nevertheless  group  theory 
effects  the  decomposition  of  the  whole  scattering  pro¬ 
cess  into  modes  of  scattering  belonging  to  the  irreduci¬ 
ble  representat i ons  of ^  .  The  dependence  on  the  row  index 
of  the  inhomogeneity  of  eq.  (9-l)  forces  us  to  solve  a  I  I 
the  systems  arising  from  the  f actor i zat i on  procedure,  in 
contrast  to  the  case  of  secular  determinants  in  which 
one  need  solve  only  one  system  per  representot i on. 

However,  as  explicitly  indicated  by  the  omission  of  the 
superscript  p,  the  T-matrjx  for  the  ^-th  mode  of  scat¬ 
tering  does  not  depend  on  p,  as  we  shall  show  presently. 
The  p- i ndcpcndcncc  o  f  &  and  J  fol  lows  from  the  very  4d<yf  i  - 

/v  A.  V.'  * 

n i t i on  of  their  matr i x  elements,  eq.  (5-7) •  In  turn  the 
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el  merits  of  the  matrices  3('CM,  %'<«),  *K«)  and  ?<V) 
are  the  symmetrized  counterparts  of  the  quantities 
and  ,  which,  as  shown  in  sect.  4b,  are  the  matrix 

elements  of  the  free  space  dyadic  Green  function  in  the 
site  and  angu I ar  momentum  representation.  The  p- indepen¬ 
dence  of  &%*)  ,5{V),  *W)  and  i s  then  a  d i rect 

consequence  of  the  inveriance  of  the  Green  function  under 
the  symmetry  operations.  Therefore  the  whole  T-matrix  for 
the  tf-th  mode  of  scatter i ng  turns  out  to  be  independent 
of  the  row  index  p.  This  circumstance  greatly  reduces 
the  computational  work  in  the  case  of  multidimensional 
repr esentat i ons .  However,  we  remark  that  the  extent  of 
the  factorization  cannot  be  stated  in  general.  In  fact 
it  depends  not  only  on  the  number  of  l- values  included 
in  the  multipolar  expansion  of  the  fields  and  on  the 
number  of  spheres  in  the  cluster,  but  also  on  the  struc¬ 
ture  of  the  symmetry  group.  The  effect  of  the  factori- 
zat  i  on  can  thus  be  illustrated  only  through  examles. 

Table  I  reports  the  order  of  the  systems  to  be  solved 

for  a  cluster  of  5  spheres  with  point  group  T  (the  CH, 

++  .  d  4 

molecule  and  the  SO^  ion  have  just  this  structure). 

Although  table  I  considers  the  case  in  which  terms  up 
to  lr=4  arc  included  in  the  expansion  of  the  scattered 
field,  from  the  discussion  of  sect.  7  we  expect  that 
terms  up  t.o  L-3  are  sufficient  to  get.  fairly  converged 
values  of  the  field.  The  useful  nos  of  group  theroy  re¬ 
quires  no  further  comment. 


Dimension  of  the  symmetrized  and  of  the  un¬ 
symmetrized  systems  for  up  to  4  for  a  cluster  of  5 
spheres  with  point  group  T^.  The  entry  U  means  unsym¬ 
metrized  while  the  other  entries  indicate  the  irredu¬ 
cible  representati ons  in  the  notation  of  ref.  25- 
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Append i \  A 

The  program  SCVSH :  a  user's  gide. 

This  appendix  is  meant  to  be  a  user's  guide  to 
SCVSH,  a  program  designed  to  generate  symmetry  adapted 
combinations  oP  Vector  Spherical  Harmonics  and  to  decom¬ 


pose  a  given  VSH  into  parts  belonging  to  the  rows  of  the 
irreducible  representations  of  a  point  group  ^  .  The  re¬ 
levant  formulas  are  reported  in  sect.  8  of  the  main  text. 
Here  we  describe  the  input  quantities  and  the  line  of 
operation  of  the  soubr out i nes . 


Al.  Input  data. 

AM  the  input  data  are  read  in  the  main  program. 
Their  sequence  and  FORMATS  are  as  follows. 

CARD  1  (3A4,SX,3I 5) 

IDSTR, IDG,N0PS,N INV,NREP 

IDSTR  1  he  name  identifying  the  structure  of  the 
cluster. 


IDG 
NOPS 
N  I  N  V 
NREP 


The  name  of  the  symmetry  group, ^  . 

Number  of  operations  in  ^  . 

Number  of  pure  rotations  in  ^  . 

Number  of  the  irreducible  represontat i ons. 


CARD  2  (1615) 

IDRCP (NP ) , ND I M(NP ) 

IDREP  Name  of  the  NP-th  irreducible  representaTji^qty 
ND I M  Dimension  of  the  NP-th  irreducible  represen- 

1 
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t  at.i  on  of 

CARD  3  (SF10.7) 

DON (NP,NR,NC0L) 

DON  Matrix  elements  of  the  first  row  of  the  NP-th 
irreducible  represent,  at  i  on .  The  index  NR  indi¬ 
cates  the  group  operation  and  NCOL  the  column. 

CARD  4  (1X,A4,5X,F3.0,3X,F4.3/3F10.7) 

ID,OM,EGA,X,Y,Z 
ID  Name  of  the  operation. 

OMfEGA  Angle  of  rotation  in  the  form  Jt^-'OM/EGA 
X,Y,Z  Direction  cosines  of  the  axis  of  rotation.  The 

program  checks  the  relation 

2  2  0 
X  +Y  +Z"=1 

CARD  5  (1615) 

LMIN, LMAX,NATOM,NFST,NLST 

LMIN#LMAX  Minimun  and  maximum  value  of  L  of  the 

VSH's  to  be  projected. 

NATM  Number  of  sites  in  the  set 

NFST, NLST  First  and  last  site  at  which  are  centered 
the  VSH's  to  be  projected 
CARD  6  (2 (3F1 0. 7 # I 5» 5X) ) 

r(i,ns),nsite(ns) 

R  x, y,  and  z  coordinates  of  the  site  NS 

NSITE  Number  oT  the  site  in  the  cluster.  In  the  main 
text  this  is  indicated  by  cA . 

CARD  7  (1615) 

NGO 
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If  NG0=0  the  program  proceeds  to  the  projec¬ 
tion  of  the  VSH's  i nc I uded  in  the  multipolar 
expansion  of  the  plane  wave.  If  NGOrO  the 
program  assumes  that  it  is  equal  to  the  num¬ 
ber  of  sites  in  the  next  set.  In  this  latter 
case  the  program  reads  a  new  set  of  cards  5 
and  6 . 


A2.  Description  of  the  subroutines. 

In  this  section  we  describe  the  subroutines 
and  functions  called  by  the  program  and  give  a  brief 
outline  of  their  mode  of  operation. 

SUBROUTINE  PRMUTE 

This  subroutine  performs  the  permutations  in¬ 
duced  on  rhe  sites  of  a  given  set  by  the  operations  of 
.  According  to  the  ma«n  text,  indeed. 


where  S  i s  an  operation  of  ^  which  for  the  present  pur¬ 
poses  is  represented  by  a  3*3  matrix.  In  the  present 
program  the  elements  of  S  are  expressed  directly  in 
terms  of  the  angle  of  rotation,  o  ,  and  of  the  direc¬ 
tion  cosines  of  the  ax  i  s,  (  ,  according  to  the  for¬ 

ma  I  as 

SuaV+(i'Xl)  =  v’X  (l-  +^U  <0 

Sic  -  X/*  (‘  -  +  ^  +(ly U1)  W  ,  =  (i'  CM  CO  )-  CO 

=  V>V*'  A~U>  ,  U*Lo)  +  \/^U>  ,  (i.  pi) 
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Th  esc  formulas  con  be  proved  by  expressing  the  matrix 


24) 


elements  of  S  by  moons  oF  the  Cay  i  ey-  K  I  e  i  r.  parameters 
as  discussed  below 

SUBROUTINE  PRJECT 

This  subroutine  uses  the  projection  operators 
of  eq.  (8-2),  to  compute  the  coefficients  of  the  symme¬ 
trized  combinations  of  magnetic  and  electric  multipoles 
included  in  the  expansion  oF  the  scattered  field.  Its 
main  ingredient  is  the  function  DL  described  below. 

The  coefficients  of  each  combination  are  stored  in  the 
(complex)  vector  COF. 

FUNCTION  DL 

This  complex  function  computes  the  elements, 
^ e  matrix  of  the  operation  S.  Instead  of 
using  the  usual  form  of  these  elements  in  terms  of  Eu¬ 
ler  angles  we  preferred  an  expression  in  terms  of 

25) 

.  Indeed,  according  to  Hammermcsh  one  has 

'Da,>  k>)  V  f  ft  -H)!(l+  H’)I(L-H’)j]^ 

MM  J*  (L-  M'-yu)J  [U'~  My*)  • 


L+M 


T"“'Vrr’r  \T  L-  b*) 


where  a.  and  b  arc  the  (£uy  I  ey-  K  I  c.  i  n  parameters  of  the 
rotation,  i.e.  the  elements  of  : 


*h.  /a-  t>  'i 
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To  compute  those  quantities  wc  recall  that  are 

the  matrix  elements  of  the  rotation  operator 

0S(„,  = 

with  respect  to  the  eigenfunctions  of  the  angu I  or  mo¬ 
mentum,  : 

/3>«}  ($)  s  <L|H'|«|»(*W'L)|L,M>s  Z*|o,*<LM'|fo.  L  |  L,M> 


Now,  it  is  well  known  that  for  L— 2  one  has 

^(><rx+/u«r3  + 

where  S'  (T«  and  O',  are  the  Pauli  spin  matrices.  Thus 

*■  1  *9  z 

w**  =  (w/.if) 


whence 


<gu=4  , 


so  that 


*(w%- 


/  .-I  v$  ^  1 1 
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s  $.  L*\  ~  ♦  i  $  A**  ~  - 


/  t  v  t 


whence 

&s  CM  ^ 

(  * 

t  -  a,(\-y)M*  ^ 


The  3x3  matrix  S  can  also  be  expressed  in  terms  of  a  and 
b  according  to  the  equation 


I  i(aVc*-  tl  -  Jl) 

\  b*t-a.  c. 


±(tx- 

{(*\t \d\£) 

Z(ac  +  Joel 


C  d  -  Cl  L 

— t(a|j+t«0 
ai/  ♦  be  i 


where  c=-b*f  d=a“r. 

SUBROUTINE  PRJPW 


This  subroutine  projects  all  the  vector  sphe¬ 
rical  harmonics  included  in  the  multipolar  expansion  of 
the  incident  plane  wave.  The  relevant  f ormu I  as  arc  cq. 
(8-2)  and  (8-6)  and  the  main  ingredient  is  still  the 
function  DL  described  above. 


t 


\:  f 


\ 


a*! titei 
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SUBROUTINE  ORTHOC 

This  subroutine  uses  the  Grahm-Schmi dt  pro- 

27) 

cedure  to  oethogona I i zc  to  each  other  the  combina¬ 
tions  generated  by  PRJCT.  A  control  is  built  in  to  re¬ 
tain  only  the  linearly  independent  combinations  which 
are  also  normalized  to  unity. 


A3*  Output  quantities. 

The  output  of  the  program  CSVSH  consists  of 
three  parts.  The  first  part  reproduces  the  names  of  the 
operations  of  and  the  corr  espond  i  ng  par  ameter  s 
described  above.  Furthermore,  the  permutations  induced 
on  the  sites  of  a  given  set  by  the  operations  of  ^ 
are  also  printed. 

The  second  part  of  the  output  gives  the  coefficients 
of  the  independent  comb i not i ons  of  magnetic  and  elec¬ 
tric  multipoles  cenctered  at  the  sites  of  the  given  set. 
The  third  part  gives  the  projection  of  the  magnetic  and 
of  the  electric  multipolcs  included  in  the  expansion  of 
the  incident  plane  wave. 

In  the  following  pages,  after  the  references,  a  list 
of  the  progam  is  included  toghet.her  with  the  input  data 
for  a  set  of  six  sites  at  the  vertices  of  an  hexagon 
(structure  of  C,H,)  with  point  group  D^.  . 
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